Inverse momentum expectation values for hydrogenic 
systems 



R. Delbourgo and D. Elliott 

School of Mathematics and Physics, University of Tasmania 
Private Bag 37, GPO, Hobart, Australia 7001 



Abstract 

By using the Fourier transforms of the general hydrogenic bound state wave func- 
tions (as ultraspherical polynomials) one may find expectation values of arbitrary 
functions of momentum p. In this manner the effect of a reciprocity perturbation 
b/p can be evaluated for all hydrogenic states. 



1 Motivation 

Many years ago and long before the elementary particle spectrum was thor- 
oughly explored experimentally, Born and Green [1] proposed the principle 
of reciprocity in an attempt to determine the mass spectrum of fermions and 
bosons. The principle flowed naturally from observed covariance of equa- 
tions of motion under the momentum-position substitution rule (P, R) — > 
(—bR, P/b), with an appropriate scale b, but it went a lot further in postulat- 
ing that the Hamiltonian was actually invariant under such a transformation. 
This was all done in a relativistic framework but it soon became apparent 
that it failed rather miserably to reproduce the mass spectrum, since all 
states were essentially harmonic overtones of a fundamental frequency; thus 
the idea was soon consigned to the dustbin of history. Recently there has 
been a resurgence in investigating the concept, not only for its elegance but 
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for its group theoretical features which incorporate the concept of a maxi- 
mum force. Much progress has been made along such lines [2] but it remains 
true that the relativistic scheme is beset with tachyonic state problems if the 
spectrum is treated along the lines of Wigner's approach of induced repre- 
sentations, ensuing from the larger 'quaplectic' group. 

Even at the classical nonrelativistic level where the idea should leave an 
imprint, the consequence for the (undamped) harmonic oscillator is that all 
frequencies are universal, which is patently absurd, and seems to indicate 
that the reciprocity concept has no future. However with the realization 
that damping with an appropriately small scale b avoids this absurdity, it 
is worth following through the idea, at least non-relativistically, for some 
familiar potentials and in particular for Coulomb-like ones about which so 
much is known. Now nonrelativistic systems placed in a 1/r potential with 
Hamiltonian 

H(P, R) = P 2 /2m - a/R (1) 

have been thoroughly studied over many years and the results are found in 
standard textbooks (see eg [3]) of classical and quantum mechanics. More 
complicated systems or small modifications, including certain relativistic cor- 
rections, can be treated by perturbation theory. If we attempt to make such 
a system reciprocity-invariant the Hamiltonian (1) is accompanied by extra 
terms b 2 R 2 /2m — ab/P. Here b is a tiny scale factor which has hardly any 
effect on atomic physics but can influence phenomena on cosmic scales [4]. 
Because b is so small the main effect of the reciprocity change lies in the b/P 
perturbation, so our aim in this paper is to determine the energy level change 
on all hydrogenic bound state levels due to it, and not just on the ground 
state as was recently done in [4]. This is the motivation for this article apart 
from its intrinsic mathematical interest. 

The neatest way to obtain the expectation value of any function of mo- 
mentum f(P) is to evaluate the Fourier transform 0(p) of the hydrogenic 
wave functions ^(x) and work out 

(f(P)) = / ^0*(P)/(P)0(P) (2) 

in the usual way. If / were purely a polynomial it would not be necessary to 
go through this process but simply work directly in coordinate space, 

</(P)) = /d 3 :rV>*(x)/HWMx), 
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because the spatial wave functions ip are very well known. However, for ex- 
pectation values of the reciprocity term where we are dealing with an inverse 
momentum, it is safest to take route (2). Nevertheless the momentum-space 
procedure (2) requires us to evaluate first the Fourier transforms, 

0(p) = J rf 3 xe- ip - x/ V(x), 

and this is technically tricky! It was first carried out in a classic paper by 
Podolsky and Pauling [5] but our treatment is slightly different and done in 
the next section. It also differs radically from Hylleraas' differential method 
[6]. The required expectation value of 1/P, the meat of the paper, is given 
in the following section; the answer is nontrivial. Interesting sum rules are 
worked out next and asymptotic estimates conclude the paper. 



2 Momentum wave functions 

We begin by quoting the bound state wave functions for Hamiltonian (1), as 
stated in the standard texts [3]: 



VWm(x) = {rd(j)\nim) = 2k 3/ \ 



(3) 

where the quantum numbers n,£,m are integers obeying n > £ > \m\ > 0, 
with k = 1/na = ma/nfi , so a connotes the Bohr radius and Y^ m are 
the spherical harmonics, orthonormal over the unit sphere. The associated 
Laguerre polynomials in (7) obey the orthogonality relations (see [7], §5.5.2) 

/ e~ t t»L»(t)L»,(t)dt = 5 nn ,T(n + v+l)/T(n + l). (4) 
Jo 

Naturally, the wave functions (3) are orthonormal over space, 



but this is difficult to demonstrate directly from (4) because of the differ- 
ent radial weight; in that respect the intrinsic n-dependence of k becomes 
significant, as shown by Dunkl [8]. 
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To be able to obtain the momentum-space wave function, 

(j) nem (\c) = (k9 k (f) k \nem) = J d 3 xe~ i]i - x ip nem (x), (5) 

we make use of the expansion of plane waves into spherical ones (see [3], eq 
(B.105)): 

e- k - = 47r^(-^(A;r)r; m (^0)r £m (^,0 fc ); j e (z) = J^-J e+1/2 (z). (6) 

By this means we arrive at (p n £m(k,9k,(pk) = ( y —i) e 'Pne(k)Ye m (9k,(fik) and the 
radial momentum wave function 



V ne (k)=4n^\ 



( " ' K : J e+1/2 (kr)(2 K rYe- Kr L 2 n % 1 _ 1 (2Kr)r 2 dr. 



^ n(n + £)\Jo V 2kr' 

(7) 

This is a formidable integral in general, except for simple cases like £ = or 

e = n-i. 

To proceed further, we make use of an integral which can be found in the 
standard texts (see [7], §3.8.3), namely 

T t+'e-fij ht) dt = ^ + Wl> + 3/2) 
and the generating function (see [7], §5.5.2) 

oo 

Y l L°(x)z v = (1 - z )- a - l e xz ^ z - l \ (9) 

i/=0 

By defining (and later on identifying v = n — £ — 1 > 0) 



= E \l 77. ■ -o , Tv777~i o, , T^s *" ( 10 ) 



{v + £ + l){v + 2£ + 1)\k 3 



we see that such an expansion in powers of z yields the momentum space 
hydrogenic functions. This new generating function is a doable integral like 

(8), 

47r 3/2 (2r ^ 



Ve(k,z) = n _^XfTk f ^e-^^J^ikr) dr 

(11) 



;i - z ) 2i + 2 V2k Jo 

8ttk (4kK) e {l- z 2 ){£ + 1)\ 

[K 2 {l + z) 2 + k 2 {l-z) 2 Y+ 2 ' 
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Since ultraspherical (Gegenbauer) polynomials arise via the generating 
series (see [7], §5.3.2), 



(1 -2xz + z 2 )- x = J2C^(x)z v ; \z\<l, A ^ 0, 



(12) 



v=0 



we are led to identify x = (k 2 — n 2 )/(k 2 + k 2 ), and derive 



oo oo 

E[^ +2 (^)-^ 2 2(^)]^ = E(^ + ^ + i)^ +1 W^7(^ + i), (is) 



v=0 



v=0 



on adopting the convention that polynomials of negative degree are identi- 
cally zero and using recurrence relations. Putting this together with (10), 
we arrive at the sought-after result for the Fourier transform, 0„£ m (k) = 
(-iyV ne (k)Y em (6 k ,(f) k ), namely 



n(n-£-l)\ (Akn)H\ 



C 



e+i 



'k 2 -K 2 



\ (n + ey. (k 2 + k 2 y+ 2 rM \k 2 + n 2 



(14) 



This will form the basis for the forthcoming calculations. [Fock's stereo- 
graphic representation [9] is also capable to reproducing (14).] 

For the present it only remains to check the normalization of these mo- 
mentum wavef unctions: 



/•OO 

/ \V nl {k)\ 
Jo 



,k 2 dk _ 32n{n-£-l)\{£\) 2 K 5 



8tt 3 



Jo 



7i(n + £)\ 
- (4k K ) 2e 

(k 2 + K 2 ) 2 ^ 



X 



a 



'k 2 ~K 2 ' 



[k 2 + K 2 



k 2 dk. (15) 



Changing the integration variable back to x — (k 2 — n 2 )/(k 2 + k 2 ), the rhs 
of (15) can be simplified to 



2n(n-£- l)\(2 e £\) 



ir(n + 



-fjl-J) e+1 ' 2 [C&_ 1 (x)]*(l-x)dx. 
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Discarding the odd term in x and using orthogonality of the Gegenbauer 
polynomials (see [7], §5.3.2), 



£ C x n {x) C x n ,{x) (1 - x 2 f~ 1 ' 2 dx = 5 nn , 



)1-2A 



vrr(n + 2A) 
(A + n)n!(r(A)) 2 ' 



we can satisfy ourselves that the rhs of (15) does indeed reduce to unity; so 
all is well for the next initiative. 



3 Momentum expectation values 

We are ready to tackle the general case, 

(f(P)U = (nim\f(P)\nim) = J \V n£ (k)\ 2 f(hk)^- 



2n(n-£-l)\(2 e £\) 



n(n + £)\ 

' (16) 

provided of course that the resulting integration is well-behaved so that f(P) 
makes sense. As we are going to be dealing with squares of ultraspherical 
polynomials, let us substitute k = KtanO, making x = — cos(2#). It then 
follows that expression (16) can be recast as 



(/( P ))^= 4W(W ^ + y !)2 / V %-2^) 2 ^(l + cos2^)[^_ 1 (cos2^)]^ 



An(n - £ - l)\(2 e e\) 2 r*/* 

(sn 

xf(hKtme)d9. (17) 
In particular, for the inverse momentum we meet the dimensionless integrals 



2n(n-£-l)\(2 e £\) 2 n 
or - 



7r(n + 



Jjl-x 2 Y[(l + x)CtU(x)] 2 dx. (18) 



We shall now show how to evaluate these for special values of £ before han- 
dling the most general angular momentum state. 
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3.1 The case I = 

To demonstrate the nontrivial nature of the problem we firstly turn to the 
spherical (S-wave) t = states. Since C T 1 t _ 1 (cos6') = sin(n6>)/ sin6> (see [7], 
§5.3.1), we must deal with 

,%k. 8 p/2 S in 2 (2n#) „ „ in 
P it Jo sin 

In order to do this (via a recursion procedure) let us define, for integer v > 0, 

K„(n)= r /2 s ™ 2{2 f cos^ede, (20) 
Jo sin 6* 



and treat the case v — first. We have 

r7r/2 

/o 



AT (n+l)-A:o(n) = /^cot #[sin 2 (2n + 2)#-sin 2 2n0] d6 = l/(2n+ 1). (21) 
•/ o 



But in the particular case n = 1 we have trivially Kq(1) — 1. It follows from 
(21) that (see [7], §1.2) 

rai 2n -i -i n -i i 

*o(n) = E = E - "2 E - = H2n+l)—Hn + l) + l (22) 

m=1 zm i m=1 m z m=1 m z z 

Another contiguity relation is 

ifi(n) -K (n) = - sin cos sin 2 2n6d9 = — -. (23) 

JO 4:Tl 2 — 1 

Combining this with (22) we obtain the final result for S-wave states, 
{(hK/P)) n0 = (4/7r)[^(n + 1/2) - 2n 2 /(4n 2 - 1) + 7 + In 4] or 

(24) 



1 8cm 

i)n0 = 



2n 2 

^(w + 1/2)- 4 ^ 2 _ x + 7 + In 4 



?'™ h 

since /i = 27r?L For large n this behaves asymptotically as 

((l/P))no ~ (8an//i)[ln(4n) + 7 - 1/2 - l/12n 2 ] + 0(n~ 3 
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3.2 The cases £ = n-l,n-2 

The case i = n — 1 is a relatively easy problem because the Gegenbauer 
polynomial collapses to unity and the integral (18) becomes quite trivial. 
We have 



S P )nn - 1 7r(2n-l)! J 1 + J 

T(n)r(n + 2) 
r(n + l/2)r(n + 3/2) 



(25) 



or 

1 _ 2na I> + !)!> + 2) 

V nn ~ 1_ r(n + l/2)r(n + 3/2) - 1 ' 

This time the asymptotic behaviour in n is 

((1/P)„„_i ~ (8na/h)[l + 3/4n + 0(l/n 2 )]. 

In a similar, but somewhat more complicated, vein we can readily treat 
the case £ = n — 2 so as to obtain 

hn _ (n + 2)r(n-l)r(n + l) _9_ /J_ 

^P /nn " 2 r(n-l/2)r(n + 3/2) 4n + U 2 

But, to proceed any further down in £, a more systematic approach is neces- 
sary. 

3.3 The general case < t < n — 1 

We can no longer avoid handling the awkward weight occurring in (18). (Had 
we been dealing with even powers of P the weight would not have been too 
troublesome, by using recurrence properties of the C^.) Realising that C^(x) 
is merely an Nth degree polynomial in x, it ought to be possible to rewrite 
it as a combination of Cj^(x) polynomials, with appropriately chosen // (to 
make the integral (18) more tractable). And indeed there is a result which 
serves this purpose (see [10], eqs 4.10.27 and 4.10.28): 

j=0 r- r(A - n)T(n + fi - j + 1) 

(27) 
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In this connection, note that for integer j > 0, lim e ^or(e + i ;')/r(e) = 5jo, so 
that (27) reduces to a triviality when X — fx. 

Let us therefore alter the weight of £ by 1/2 up and down by expressing 



[(n-/-l)/2] 
3=0 



■e+i/2 

n-£-l-2j 



Vn-t-l)/2] 

(*)= E T^-fi^W, (28) 

J"=0 



where 



(n - 2j - 1/2) T{£ + 1/2) T(j + 1/2)I> - j) 

Pjnt = 



Ijni 



j\ T(£+l) r(l/2)r(n-j + l/2)- 

(n - 2j + 1/2) T(£ + 3/2) T(j - 1/2)I> - j) 



T(£+l) r(-l/2)r(n-j + 3/2)- 



(29) 
(30) 



Then, using the orthogonality property of the ultraspherical polynomials, we 
end up with the single summations: 



r l [(n-<-l)/2]/ cy-ln \ 

/ji-.V^Ml^ g (^) 



(ii-2j-l/2)(ii-£-2j-l)! 



and 

-l 



[(ra ^' 1)/2 V r%u V 7rr(n+£-2j) 
J ' 



[(n "'- 1)/21 ^- 1 -SnA 2 7rr(n+£-2 J + 2) 



^ V r (^+ 3 / 2 )/ (n-2j + l/2)(n-£-2j-l)! 
Applying these to (18) we obtain 



5« 2n(n-l-l)! / T(£ + l) V [(ra ^ 1)/2] r(n + £ - 2j) 

{ P )n£ ~ (n + £)\ [r(£ + l/2)) ^ V{n-£-2j) 



x 



2(3 2 



jnt 



( n + e _ 2 j)(n + £ + l-2j) 1 ] ne 



n-2j- 1/2 



(2£+l) 2 (n-2j + l/2) 



(31) 



Substituting the expressions for f3 and 7 from (29) and (30) respectively, we 
can finally reduce the answer to 

Tin 2n{n-£-l)\ [{n ^ )/2] f T{j + l/2)T{n-j) \ 2 F(n+£-2j) 

{ T U ~ n(n + £)\ ^ 



x 



(271-1-4?)- 



J=0 K T(j + l)T(n-j + l/2)J V{n-£-2j) 
(n+£-2j)(n+l/2-2j)(n+£+l-2j) 



(2j-iy(2n-2j + lf 



(32) 



£\n 


1 2 3 4 5 6 





32 256 2144 1024 85088 1172224 
3 15 105 45 3465 45045 


1 


128 256 512 57088 809344 
15 21 35 3465 45045 


2 


4096 16384 299088 21856256 

525 m m — wffi- 


3 




4 


2205 b& 

72765 


5 


297297 



Table 1: Calculated values of {2nhK,/P) n £ for principal quantum number n 
running from 1 to 6, as given by (32). 

We have not succeeded in simplifying this any further, except for the earlier 
simple cases. In Table 1, we have provided the first few expectation values of 
2tiTik/ P (for n and I up to 6) as derived from (32). Remember that k = 1/na, 
where a is the Bohr radius. 



4 Sum rules and integral representation 

There exists an interesting set of sum rules for {hK,/P) n £ which hail from 
a particular form of the addition theorem for Gegenbauer polynomials (see 
[7], §5.3), and which produce a neat integral representation. Start with the 
particular 4-dimensional case 

C;Ucos 2 + sin 2 0cos^) = g ~ "L^f 7TT T ~ ( 2 sin ^ x 

JT Q 1(71 + 1+1) 

lCi + _U(cos6)] 2 Pe(cos^). 

On putting cos 6 = x, we find that the rhs has a close connection with the 
rhs of eq (18) and so recognize that 



/ (1 + xfC^x 2 + (1 - x 2 ) cos^) dx = — ]T(2£ + l)P e (cos^)(hK/P) rd . 
J-i In e=Q 



(33) 



10 



fj l~ On r 1 

J2(2£+l)P e (y)(—) ne = — (1 + x 2 )U n ^(x 2 + (1 - x 2 )y) dx. (34) 
tr n r 7r J-i 



Since C^_ 1 (y) = U n -i(y), we can simplify this result to 

g^. , .^k. 2n /-i 

To obtain the sum rules, simply set y — 1 and y — — 1: 

^/ „ 2n 2 f 1 . , 16n 2 . , 

£2*+l )(—) nl =— (l + x 2)dx = —, (35) 
S P vr 7-i 3tt 

£(2€ + l)(-l)<<— )nt=— (l + x 2 )U n ^(2x 2 -l)dx. (36) 

In order to evaluate the rhs of eq (36) first define J n = j\ U n {2x 2 — 1) dx. The 
rhs of (36) can then be rewritten as (J n +6J n _i + J n _ 2 )/4, so it only remains to 
determine J n . This can be done via an easily established recurrence relation, 
namely J n + J n _i = 2/(2n + 1). Hence we deduce that 

2J n = ijj(n/2 + 5/4) - ijj(n/2 + 3/4) + (-l) n vr. 

Therefore 



,Hk, n 



xi(2£+i)(-i/Q n£ =- 



(37) 

Finally, on multiplying (34) by a Legendre polynomial in y and inte- 
grating over (-1,1), we obtain a neat double integral representation for the 
expectation value. Thus 



, tin . n r 1 2 . 



j\ e (y)U n ^(x 2 + (l-x 2 )y)dy 



dx. (38) 



One may verify in particular cases that this produces the same results as the 
series (32) but, more significantly, it allows us to obtain a series representation 
of (hK,/P) n £i which is different from (32). 

To arrive at this new result, we make use of the fact (see [7], §5.7.2) that 



^1(n-j-l)!W(j|3/2)' 
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and notice the factorization, 1 — (x 2 + (1 — x 2 )y) = (l—x 2 )(l—y). Thus (38) 
factorizes into two parts. The integral over x is easily done, but the integral 
over y is a bit harder (requiring Rodrigues' formula and an integration by 
parts). Carrying out the necessary manoeuvres we end up with 

Hk (-I)' n(£ + j + 2)(n + £ + j)\[(£ + j)!] 2 

[ p )n£ p (n-e-j-i)\(2e + j + i)\j\r(e + j + s/2)r(e + j + 5/2) ■ 

(39) 

This is more compact than the series (32). 



5 Asymptotic behaviour 

There are 3 regimes to consider for large n which we shall look at in turn. 
5.1 l/n«l 

Returning to eq (18) we may replace F(n — £)/Y{n + £ + 1) asymptotically 
by rT 2t ~ x and use the approximation (see [7], §5.3.3), 

, „„x n e cos(2n# - ir(£ + l)/2) , n n 

This tells us that 

,hK. , 16 W2- ec os 4 ^cos 2 (2^-7r(£+l)/2) ,„ 
\-5-)nt ~ lim — / -^r Q d9, 

P e^0+ 7T Je Sin 26* 

assuming the integral exists. Hence we obtain the asymptotic difference 
between two £ values differing by two: 

tin. ,Tik 16 fx/i-t cos 4 6 
( — ) n e+i — \-^r)ne-i ~ hm — / dO -x 

\p/nl+L \ pint e^Q+ n J e s i n 2# 



cos 2 (2n9 - n(£ + 2)/2) - cos 2 (2n# - n£/2) 



0. 



But we have already shown in (24) that (Hk/ 'P) n i ~ 4^(n+l/2)/7r for £ = 0, 
from which we conclude that 

(^>„ £ ~4^(n+l/2)/7r (40) 

for n ^> 1 and modest values of £. This result diverges logarithmically like 
(41ogn)/7r as n — > oo. 
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5.2 £/n near 1 

In this case return to the series (32) and put n = I + 1 + 5 where 5/n <ti 1. 
Making use of the asymptotic approximation, 



T(z + a) a _ b 
hm t4- ~ z 



z->co Y(z + b) 
the sum collapses to 



l + ^-(a-b)(a + b+l) + 0(\) 



(41) 



TlK 



-1-8 



1 



7T 



[5/2] 



E(2n) 



-2J. 



r(i+<s) /T0'+i/2) > 



2r 



i=o 



r(i+5-2j)Vr(j+i) 



(2j-i) 



1 + *A 



n 



Thus j = is the dominant term in the expansion; this is multiplied by a 
subdominant factor of order 1/n via (32), producing the final estimate, 



Tin 

\~p)nn-l-S ~ 



An n 2 



(42) 



This agrees with the results of section 3.2 where the cases 5 = and 1 were 
studied. 



5.3 £/n is finite 

This regime is the trickiest to deal with, interpolating between (40) and 
(42) as it does. Although we have not succeeded in obtaining the analytical 
dependence on A = £/(n — 1), one may readily establish numerically that 
the sum (39) tends to a constant in the limit of large n and increases as A 
approaches zero; for instance when A = 1/2 the value is 1.975; when A = 1/4, 
the value is 2.88; when A = 1/8, the value is 3.77, etc. This dependence would 
appear to be logarithmic as is indicated by (24). 



6 Conclusions 

In this paper we have principally concentrated on the evaluation of (1/P) 
- a challenging problem — because it has special significance for Born reci- 
procity, but the methods we have used can no doubt be extended to general 
functions of momentum via eqs (16) and (17). In fact such calculations have 



13 



been carried out in [11] for expectation values of P N and logP to which the 
reader is referred. At any rate, the conclusion of our analysis is that the 
1/p perturbation has an increasingly disruptive effect on the higher n states 
having the lowest £ (the dependence is logarithmic in n). This hails from 
the full effective potential [4] which shows a maximum below zero energy, at 
E = bL — 2a\^bL, where L signifies the classical angular momentum. 
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